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Abstract 

The formalism for the neutrinoless (^~,e + ) conversion is investigated in detail and 
the relevant nuclear matrix elements for light intermediate neutrinos in the case of 
27 Al(/j,~, e + ) 27 Na are calculated. The nucleus 27 A1 is going to be used as a stopping 
target in the MECO experiment at Brookhaven, one of the most sensitive probes 
expected to reach a sensitivity in the branching ratio of the order 1CP 16 within the 
next few years. The relevant transition operators are constructed utilizing a variety 
of mechanisms present in current gauge theories, with emphasis on the intermediate 
neutrinos, both light and heavy, and heavy SUSY particles. The nuclear wave 
functions, both for the initial state and all excited final states are obtained in the 
framework of Is — Od shell model employing the well-known and tested Wildenthal 
realistic interaction. In the case of the light intermediate neutrinos the transition 
rates to all excited final states up to 25 MeV in energy are calculated. We find that 
the imaginary part of the amplitude is dominant. The total rate is calculated by 
summing over all these partial transition strengths. We also find that the rate due 
to the real part of the amplitude is much smaller than the corresponding quantity 
found previously by the closure approximation. 



PACS number(s): 23.40.Bw, 23.40.-s, 14.60.Pq, 21.60.Cs 







1 Introduction 



Modern gauge theories (grand unified theories, supersymmetry, etc.), which go beyond the 
standard model (SM), predict a great number of processes which violate the lepton and/or 
lepton family (flavor) quantum numbers IQ-0. One can distinguish three categories of 
such processes: the purely leptonic (// — > e + 7, \i — > 3e, etc.), the semileptonic hadron 
decays [K\ — > /i ± e ± , etc.), and the semileptonic which take place in the presence of nuclei. 
Among the most interesting such processes are those which take place in a muonic atom 



One exotic such process is the muon-to-positron conversion, 



H~ + (A, Z) ^e + + (A,Z-2) 



which violates the muonic (L^), electronic (L e ) and total lepton (L) quantum numbers 
Another exotic process is the muon-to-electron conversion, 



17-27 



(A,Z)^e- + (A,Z) 



(2) 



which violates only the lepton-family quantum numbers L M and L e 

Both of these processes can experimentally be studied simultaneously, since both of 
them have the same intrinsic background and the same initial state (a muon at rest in the 
innermost IS orbit of a muonic atom). In the present work we will focus our attention on 
reaction ([]]). 

In recent years, continuous experimental efforts have been devoted for the measure- 
ment of the branching ratio R fie + defined as the ratio of the (fi~,e + ) conversion rate 
divided by the total rate of the ordinary muon capture reaction pOL i.e. 



R 



i>- - e + )/r(^- 



(3) 



Up to now only upper limits have been set and the best limit found is for the 48 Ti nucleus 
at TRIUMF and PSI O, O ITS, 



R 



yielding the values 
,-12 



< 4.6 x 10" 







and 



R»e+ < 4.4 x 1(T 12 H 



respectively. This limit could be further improved by future experiments, at PSI (SIN- 
DRUM II experiment), which aims to push the sensitivity of the branching ratio R^ e + to 
10~ 14 , and at Brookhaven (MECO experiment) [28[] with expected sensitivity about three 
to four orders of magnitude below the existing experimental limits H28|, |2l 

Traditionally the exotic /i 
(like 48 Ti and 63 Cu 



e processes were searched by employing medium heavy 



|]~3|, 16] or very heavy (like 208 Pb and 197 Au) |14, 15, 16] targets. For 
technical reasons the MECO target has been chosen to be the light nucleus 27 A1 f28| . 
The MECO experiment, which is planned to start soon at the Alternating Gradient 
Synchrotron (AGS), is going to employ a new very intense \T~ beam and a new detector 
|28| . The basic feature of this experiment is the use of a pulsed \i~ beam to significantly 
reduce the prompt background from 7r~ and e~ contaminations. 

The best upper limit for the yT — > e~ conversion branching ratio R^e- set up to the 
present has been extracted at PSI (SINDRUM II experiment) as 



R*e- < 5 -° x 10 ' 



13 



for 19 'Au target [16] . 



(4) 
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For the 48 Ti target the determined best limit is 



R 




< 6.1 x 10 



13 



[13], 



(5) 



while for the 208 Pb target the extracted limit is 



R 



,Pb 

7 te 



< 4.6 x 10 



ii 



[14]. 



(6) 



Processes (1) and (2) are very good examples of the interplay between nuclear and 
particle physics in the area of physics beyond the standard model. Moreover, the e + ) 
conversion has many similarities with the neutrinoless double /3-decay (Oz//3/3), 



and especially with its sister electron to positron conversion || , which violate the lepton- 
flavor (L e ) and total lepton (L) quantum numbers. Both reactions (P and (0) involve 
a change of charge by two units and thus they cannot occur in the same nucleon. Both 
of them are forbidden, if lepton number is absolutely conserved. One can show that, if 
either of these processes is observed, the neutrinos must be massive Majorana particles. 
In spite of the many similarities, however, these double charge exchange processes have 
some significant differences which can be briefly summarized as follows: 

(i) Due to the nuclear masses involved, neutrinoless double beta decay can occur only in 
specific nuclear systems for which single beta decay is absolutely forbidden, due to energy 
conservation, or greatly hindered due to angular momentum mismatch. These systems, 
with the possible exception of 48 Ca, have complicated nuclear structure. Furthermore, 
the neutrinoless double beta decay can lead mainly to the ground state and, in some cases, 
to few low-lying excited states of the final nucleus. Such constraints are not imposed on 
process (|1|), due to the rest energy of the disappearing muon. 

(ii) From the corresponding experiments, in conjunction with appropriate nuclear ma- 
trix elements as input, one may extract lepton violating parameters, which depend on 
flavor. Thus, in the framework of the neutrino mixing models, the amplitudes of neutri- 
noless (3/3 decay and e + ), for leptonic currents of the same chirality, are proportional 
to some combination of the neutrino masses, which are different from each other. The 
same is true in the case of the mass independent lepton violating parameters r] and A 
appearing in the amplitude, if the leptonic currents are of opposite chirality. One does 
not know a priori which flavor combination is favored. 

(iii) The long wavelength approximation does not hold in the case of (fi~, e + ) conver- 
sion, since the momentum of the outgoing e + is quite high. Thus, the effective two-body 
operator responsible for the (/i~,e + ) conversion is strongly energy dependent and more 
complicated than the corresponding one for the 0^/3/3 decay. On the other hand, one can 
in this case choose a target, consistent with the standard experimental requirements, so 
that the nuclear structure required is the simplest possible. 

(iv) Neutrinoless double beta decay has the experimental advantage that there are 
no other fast competent channels. The only other channel is the two-neutrino double- 
beta decay, which, however, is also of second order in weak interaction and kinematically 
suppressed. 

In any case we view the two processes, (/i~,e + ) conversion and OvjSP decay, as pro- 
viding useful complementary information. As is well known in the case of neutrinos one 
needs the following information: a) the mixing matrix, b) the three eigenmasses and c) 



(A,Z) -> e" + e~ + (A, Z + 2) 



(7) 
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the two independent Majorana phases. These phases are present even in a CP conserving 
theory since they can take values and tt. Study of CP violation in the leptonic sector is 
beyond the goals of the present paper. The interested reader is referred to the literature, 
e.g. see j31 . 



The neutrino oscillation data can, in principle, determine the mixing matrix and the 
two independent mass-squared differences (e.g. Sm^ = m 2 . — mf, Sm^ = m\ — mf). This 
has been done in a number of papers, see e.g. the recent work by Haug et al. |32| and 
references therein. They cannot determine, however, the scale of the masses, e.g. the 
lowest eigenvalue mi and the two relative Majorana phases. Once mi is known, one can 
find 7722 = [^tti^ -l-m 2 ] 1 / 2 and 7773 = [Sm^ + 777 2 ] 1 / 2 . The masses can be chosen positive by 
a proper redefinition of the neutrino fields, i.e. by absorbing the sign into the Majorana 
phases. If the mixing matrix is known, the mass mi can be determined from triton decay 
as follows: 

(m u ) ee = m v = | £ KP^ 1 ' 2 > (8) 

3=1 

since the Majorana phases do not appear (this experiment cannot differentiate between 
Dirac and Majorana neutrinos). In the above expressions U is the mixing matrix. 

The two lepton violating processes provide two independent linear combinations of the 
masses and the Majorana phases. In fact 

3 

{m u ) ee = (m v ) = | ^2 U ej U ej e' Xj rrijl (9) 

3=1 

and 

KW = IEW^"» J |. ( 10 ) 

3=1 

It is clear from the above expressions that only the two relative CP phases are measurable. 

So, the three experiments together can specify all parameters not settled by the neu- 
trino oscillation experiments. So, if possible, all three of them should be pursued. 

Strictly speaking (/j~,e + ) and neutrinoless double beta decay should be treated as 
two-step processes by explicitly constructing the states of the intermediate (A, Z ± 1) 
systems. It has been found [^3|, however, that, for neutrinoless double beta decay, since 
the energy denominators are dominated by the momentum of the virtual neutrino, closure 
approximation with some average energy denominator works very well. We expect this 
approximation to also hold in the case of (/i~,e + ) conversion to sufficient accuracy. We 
will, therefore, replace the intermediate nuclear energies by some suitable average. By 
summing over the partial rates of all allowed final states of the nucleus (A, Z — 2) we 
obtain the total rate. This will then be compared to that obtained by invoking closure 
T9fl with some appropriate mean energy (Ef) of the final states. 



So far, theoretically the (/i , e + ) process has been investigated [T9|, [24j] on the exclusive 
reactions 

40 Ca + /i- ^e + + m Ar(gs) (11) 

58 Ni + fT -> e + + 58 Fe(gs) . (12) 

In these studies the partial g.s. — > g.s. transition rate was calculated by performing 
microscopic calculations of these nuclear matrix elements. On the other hand, the total 



3 



transition strength to all final states (inclusive process) was estimated along the lines of 
closure approximation and ignoring 4-body terms [[HJ . 

In the present article we apply the shell-model techniques in the investigation of the 
(/i~,e + ) conversion reaction 

27 Al + fi~ -> e+ + 27 Na . (13) 

To this end, we construct all the needed wave functions and calculate the rate not only 
to the ground state but also to the excited states of the final nucleus. As a first step we 
limited ourselves to the mass mechanism for light neutrinos. So, only the multipolarities 
associated with the Fermi and Gamow- Teller type operators were relevant. For practical 
reasons we had to limit ourselves to states lying below some excitation energy (« 25 
MeV) of the final nucleus 27 Na. The distribution of the strengths for the most important 
multipolarities versus the excitation energy of the daughter nucleus in the reaction (p~3|) 
is also studied. 

The paper is organized as follows. In sect. 2, an extensive presentation of the relevant 
expressions occurring in the formal description of the jjT —>■ e + transition operators is 
given. In sect. 3, we deal with the expressions of the branching ratios. In sect. 4, we 
discuss the evaluation of the inclusive /i~ — > e + matrix elements by means of explicit 
construction of the needed nuclear wave functions in the framework of the s — d shell 
model. In sect. 5, the results obtained for the transition matrix elements in the case of 
27 Al(fi~, e + ) 27 Na are presented and discussed. Also the spreading of the contributions due 
to the occurrence of various multipoles are described. Our conclusions are summarized in 
sect. 6. 

2 Brief theoretical formulation of the yi~ — » e + con- 
version operators 

2.1 Effective ji~ —> e + conversion Lagrangian in Gauge models 

^From the particle physics point of view, processes like \i~ — > e ± conversions are forbidden 
in the SM by total- lepton and/or lepton-flavor (muonic and electronic) quantum number 
conservation. They have been recognized long time ago as important probes for studying 
the lepton and lepton-flavor changing charged- current interactions which go beyond 

the standard model. 

There are several possible elementary particle mechanisms which can mediate the 
lepton- violating process (Q). The mechanisms which have been studied theoretically are: 

(i) Those mediated by a massive Majorana neutrino. In this case we have two possi- 
bilities. 

1) The chiralities of the two leptonic currents are the same. Then, the amplitude 
in the case of light neutrinos is proportional to some average neutrino mass or to some 
average of the inverse of the neutrino mass, if it is heavy. 

2) The chiralities of the leptonic currents are opposite. Then, the amplitude is not 
explicitly dependent on the neutrino mass. It vanishes, however, if the neutrinos are not 
Majorana particles. This mechanism is significant only in the case of light neutrinos. 

(ii) Those accompanied by massless or light physical Higgs particles (majorons). 

(iii) Those involving more exotic intermediate Higgs particles. 
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(iv) Those mediated by intermediate supersymmetric (SUSY) particles. 

From a nuclear physics point of view one has to be a bit careful when the intermediate 
particles are very heavy. The elementary amplitude is constructed at the quark level, 
but the calculation is performed at the nucleon level. If, in going from the quark to the 
nucleon level, the nucleons are treated as point-like particles, the nuclear matrix elements 
are suppressed due to the presence of short range correlations. To avoid this suppression, 
a cure has been proposed || |34j which treats the nucleons as composite particles described 



by a suitable form factor. A different approach is to consider mechanisms, which involve 
particles other than nucleons in the nuclear soup. Such are, e.g., mechanisms whereby 
the processes ([!]) and (0) are mediated by the decay of the doubly-charged A ++ (3/2, 3/2) 



resonance 0, |23|, [27| present in the nuclear medium, i.e. 

\f + A ++ -> n + e + (14) 

This process, however, to leading order does not contribute to + — > + (3(3 decays, but it 
may contribute to (/i~,e + ). The other process is induced by pions in flight between the 
two nucleons [||, [36[ according to the following modes: 

(a) The 1-pion mode represented by the reactions: 

^ + p — * n + t\~ + e + , n~ + p — >n (15) 

where the protons and neutrons are bound in the nucleus. 

(b) The 2-pion mode represented by the reactions 

P — > n + 7T + , 7T + /i~ — > 7T~ + e + , 7T~+£>— > 11 (16) 

2.2 The transition operators at nuclear level 

The gauge models mentioned in the previous subsection give rise to a plethora of effective 
transition operators. Their essential isospin, spin and radial structure is given as follows. 

1) The isospin structure is quite simple, i.e. of the form r_(z)r_(j) where i and j label 
the nucleons participating in the process. 

2) The spin structure is given in terms of the operators 

W S i{ij) = 1 (Fermi) (17) 

Ws2(ij) = °i ■ °~j (Gamow — Teller) (18) 

Wsa{ij) = 3(<Tj • f ij )(a j ■ f iS ) - a { ■ crj (Tensor) (19) 

W A i{ij) = %<Ji x Oj (20) 

W A2 (ij) = <Ti - Gj (21) 

The operator is determined below. 

3) The orbital part can be expressed in terms of the quantities: 

(a) The momentum of the emitted positron (p e ) obtained from the kinematics of the 
reaction (fj). One finds that 

Pe = \p e \ = rn^ - t b + Q - E x (22) 
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where Q = M(A, Z) — M(A, Z — 2) is the atomic mass difference between the initial, 
(A, Z) and final, (A, Z — 2) nucleus, e& is the binding energy of the muon at the muonic 
atom (eft ~ 0.5 MeV for 27 A1), E x is the excitation energy (E x — Ef — E gs ) of the final 
nucleus and m M is the muon mass (m M = 105.6 MeV). 

(b) The relative (r^) and center of mass (Ry) coordinates which are written as 

and ^ 

I I 

The radial part of the operator contains: 

(i) the spherical Bessel functions ji(p e fij/2) and jciPeRij) resulting from the decom- 
position of the outgoing positron and 

(ii) a function f(r) of the relative coordinate (r = ry) given by: 

f(r) = ^F(r)V cor (r), (23) 
r 

where the constant Rq represents the nuclear radius (Rq = \.2A 1 ^). The function \l/ cor (r) 
is some reasonable two-nucleon correlation function, e.g. of the type [[| 

V cor (r) = l-e- ar \l-br 2 ) (24) 

with a = \.\jmT 2 and b = 0.68/m -2 . 

We mention that, strictly speaking in the above radial function the muon wave function 
must also be included. In fact, if muon is considered as relativistic particle an additional 
lepton-spin dependence appears in the transition operator. In the case of the light nucleus 
27 Al studied in the present work, however, the muon is in the Is atomic orbit. Its wave 
function varies very slowly inside the nucleus and thus it can be replaced by its average 
value. Anyway this average value drops out if the same approximation is assumed for the 
/i~-capture, i.e. for the denominator of the branching ratio R^ e +- 

As we have already indicated, the radial function F(r) depends on the specific mech- 
anism assumed for the /i~ — * e + conversion process to occur. The following cases are of 
interest . 

(i) In the case of light Majorana neutrinos, when the leptonic currents are left-handed, 
F(r) takes the form 



F(r) = - / dx + - / dx (25) 



7T Jo x — Ot + 16 71 Jo X + 5 

The quantities 6 e and a are given in terms of the nuclear masses and the average excitation 
energy of the intermediate states, (E xn ), as 

5 e = [(Exn) + M(A, Z-l)- M(A, Z) + p e ]r 

a = K + M(A, Z) - M(A, Z-l)- (E xn )}r 
Note that 5 e depends on the positron momentum. The first term of F{r) in Eq. (EH) can 



be written as 

2 r°° sinx , 2 f°° sinx , . 

-dx = — F I dx — i2sma (26) 



2 r c 

7T JO 



x — a + te 7r Jo x — a 
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As can be seen, the amplitude has now an imaginary part, a fact that was missed in 
earlier calculations The principal value integral can be written in an equivalent form 
as follows: 

2 f°° sinx , 2 f°° sinx , . 
— r I ax = zcosa — 1 H — a / — -ax {27) 



7r Jo x — ol tt Jo x{x + a 

The latter expression is more convenient for numerical integration techniques. By com- 



bining Eqs. (^j)-(27) we can write F(r) with separate real and imaginary parts as 

2 f°° r a 1 1 , / 
2cosa — 1 H — / sinx — H ax — i2sina . (2£ 

7i Jo l x(x + a) x + o P j 



F(r) 

71 JO 



It is worth remarking that, in the case of 0i>(3(3 decay a ~ therefore F(r) = 1. This 
simplifies quite well the calculations in the 0^(3 (3 decay process. 

(ii) In the case of light Majorana neutrinos, when the leptonic currents are of opposite 
chirality we have F(r) — > F (r) = r(d/dr)F(r). The same situation occurs in the context 
of R-parity violating supersymmetric interactions mediated by light Majorana neutrinos 
in addition to other SUSY particles. 

(iii) For heavy intermediate particles, e.g. heavy Majorana neutrinos, we will examine 
two modes: 

1) Only nucleons are present in the nucleus. Then the function F(r) reads 

F{r) = ^-^-x A (x 2 A + 3x A + 3)e- XA , x A = m A r (29) 
48 m e m p 

with m e , m p the masses of electron and proton respectively. It should be mentioned that 
the above expression was obtained for neutrino masses much heavier than the proton 
mass provided that the nucleon is not point like. The particular expression holds, if the 
nucleon is assumed to have a finite size adequately described [|^] by a dipole shape form 
factor with characteristic mass m A taking the value m A = 0.85GeV/c 2 

2) The process is mediated by pions in flight between the two interacting nucleons. 



Then one distinguishes two possibilities [y, [39 



(a) The 1-pion mode Eq. (|T^). In this mode F(r) is replaced by F[ w , % = GT (Gamow- 
Teller), T (Tensor) where 

F<f (x) = a ln e~ x , F^(x) = a ln (x 2 + 3x + 3)e~ x /x 2 (30) 

with x = x-n = m^r (m^ denotes the pion mass) and = 1.4 x 1CT 2 . In this case the 
radial functions are the same with those entering the neutrinoless double beta decay. 

(b) The 2-pion mode Eq. (0). Now the radial functions are very different from those 
appearing in 0z/ (3/3 decay, since the momentum carried away by the outgoing lepton is 
not negligible. It is, however, obtained from those entering the neutrinoless double beta 
decay, via the substitution: 

F(rH(x e /2) - f'm ~ l/2)x e )FUm ~ 04 + 4] 1/2 R (31) 
Jo 

(x e = m e r) where F^, i = GT, T are given by [|J 

Fg T (x) = a 27T (x-2)e~ x , ifjx) = a 2n (x + l)e~ x (32) 
with a 2 n = 2.0 x 10~ 2 . 
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2.3 Irreducible Tensor Operators 

In this section we are going to exhibit the structure of the various irreducible tensor 
operators relevant to our calculations. We characterize them by the set of quantum 
numbers /, £, A, /i, A, L, S, J, some of which may be redundant in some special cases. Thus 
I and C refer to the multipolarity of the outgoing lepton in the relative and center of mass 
systems. A is the orbital rank of the operator in the relative coordinates (the corresponding 
rank in the CM system is /x). L is the total orbital rank, S is the spin rank and J the total 
angular momentum rank. Finally A is the rank of the spherical harmonic describing the 
momentum of the outgoing lepton. The latter couples to zero or 1 with the J-rank, i.e. 
(A, J)k, k — 0, 1. Some details on how these operators are combined to give the nuclear 
matrix elements will be discussed in the Appendix. 

We will begin with operators appearing when the chiralities of the two leptonic currents 
involved are the same. This covers the case of minimal left-handed extensions of the 
standard model. 

One encounters Fermi- type operators, ftp, of the form 



i<j T - ^) T - U ) / ( r U )jl ( ^ )3C (PeRij ) 



A 



l,fi = £, 



5 = 0, J = L 



(33) 



(A is redundant). The Gamow- Teller operators, ^gt, are similarly written as 



GT 



Ei<j T - 



4ttF A (fij ) <g> V4^Y» (R i:j ) <g> ( - y/3) fa a,] 



A = /, fi = C, S = 0, J 



(A is redundant). Note that 



<7j • Gj = -a/3 [<7j 

The first spin antisymmetric operator is 



a 



no 



(34) 
(35) 



Q A1 — 



i)r-{j)f( r ij)ji( 
A = I, ji = C, S = 



E«<j 
4irY x (r. 



<g> (->/2) [<Ti <g> <Jj 

1, J = L,\L±1\ 



(36) 



(A is redundant). Note that 



[-V2) [ji 



(37) 



The second spin antisymmetric operator is 



^A2 = 



47ry A (r^) <g> V^Y^R^l <g> (^ - <tj) 
A = /,/i = £, 5 = 1, J = L,\L±1 



(38) 
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(A is redundant). Note that each operator must be overall symmetric with respect to 
interchange of the particle indices. So, in those cases in which the spin operator is of rank 
unity, I must be odd. In the special case of + — > + neutrinoless double beta decay, only 
the Fermi and Gamow- Teller operators occur. 

We are now going to consider the case in which the theory contains both R (Right) 
and L (Left) currents. If both currents are right-handed the above results hold, but 
the relevant neutrinos are heavy. We thus need only consider the case in which we 
have L — R interference in the leptonic sector. This may be important in the case of 
light neutrinos. As we have already mentioned, this also occurs in the context of R- 
parity violating supersymmetric interactions, which, in addition to other SUSY particles, 
involve intermediate light Majorana neutrinos. The amplitude now is proportional to the 
4-momentum of the intermediate neutrino. The time component has a structure similar to 
the above with a different energy dependence. The corresponding operators are indicated 
by putting a "prime" over the corresponding ones for the mass term. This point will 
not be further pursued, since their contribution is suppressed. Its space component, after 
the Fourier transform, gives an amplitude proportional to the gradient of the Fourier 
transform of the previous case. We thus get the above operators, to be denoted by fl F 
Q GT , and Q A2 (associated with the term linear in the spin), with F(r) replaced by F (r). 
Now the overall operator has an overall rank of a vector, obtained by the coupling of the 
two operator ranks J and A. In this case, in addition to operators of the above form, we 
encounter an operator of spin rank two, which is of the form: 

ti T = T-{i)r-.{j)f{r ij )3i{^-)3cipJtii) 
^Y x (f tj ) ® y/birY^Rij)] ® [<Ti ® (Tjf 
\= \l±l\,n = C, S = 2, J = L,\L±1\,\L±2\ (39) 

(A is redundant). 

As it has already been mentioned, in the case of heavy intermediate particles one may 
have to consider pions in flight between nucleons. Then one encounters only Gamow- Teller 
and tensor operators except that now the radial part is different (see Eq. (j3Cf)-(l3~2|)). 

In the special case of + — > + neutrinoless double beta decay mediated by light 
neutrinos one can invoke the long wavelength approximation. Thus, to leading order one 
finds (up to normalization constants and possibly factors of p e ) the familiar operators: 

^ = E^W^(J')/N (Ferrrn) (40) 

^gt = T -(i) T -U)f( r ij) a i ' °i {Gamow - Teller) (41) 
ti A 2 = E T-WT-Wf'fa) (o-i - a^.iifij x Rij) (42) 
^'t = T -( i ) T -ti)f ( r ij) [ 3 ( a i ' ' ~ a i ' a j) (Tensor) (43) 

i+3 
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3 Branching ratio 



The branching ratio R^ e + of the (p~,e + ) reaction defined in Eq. (|3|) contains the LFV- 
parameters of the specific gauge model assumed. These parameters, are entered in i? Me + 
via a single lepton- violating parameter n e ff. Under some reasonable assumptions these 
parameters can be separated from the nuclear physics aspects of the problem. As has 
been pointed out p4j, the branching ratio -R Me + takes the form 



with T] e ff = (m v ) / m e The parameters p and r] e ff, which depend on the gauge model 
adopted and the mechanism prevailing, are expected to be very small due to the fact 
that \i~ — > e + conversion is a lepton violating second order weak process 0. In this 
definition, the total muon capture rate has been written in terms of the well-known 
Primakoff function Jpr(A, Z) ||0| which takes into account the effect of the nucleon- 



nucleon correlations on the total muon capture rate. |A^i^/| 2 denotes the square of the 
partial transition nuclear matrix element between an initial \i) and a final |/) state. This 
can be written as 

|Xwl 3 =9^rr £ \iJfMf\n\JiMM 2 (45) 

ZJi 1 M f Mi 

In our case \i) = \g.s.), i.e. the ground state of the initial nucleus. The summation in 
Eq. P^D runs over all states of the final nucleus lying up to ~ 25 MeV. (|/) = |JJ). We 
consider the final nuclear states as having well-defined spin (J) and parity (n). The index 
p counts the multipole states). 

For the Fermi and Gamow- Teller contributions, which are expected to be the most 
important contributions, the square of the matrix element |.Mj^/| 2 is written as 



\Mi 



*f\ 



|2 



2Ji 




i{gs)) - (f\\Q GT \\i(gs)) 



(46) 



where fy and Ja are the usual vector and axial vector coupling constants (/a//v — 1-25). 
By combining Eqs. ( fPj) and (|46|) we see that, for the evaluation of the branching ratio 
R^ e +, we have to calculate the reduced matrix elements and (/H^Grll*) for 

\i) = \g.s.) and |/) any accessible state of the final nucleus. In the present work these 
states have been constructed (for 27 A1 and 27 Na systems) in the framework of the shell 
model as is described in the next section. 



4 The shell model nuclear wave functions 

The evaluation of the reduced matrix elements (/| \flp\ \i) and (/| | fieri \i) f° r the p~ — > e + 
conversion requires reliable nuclear wave functions. These were obtained in the framework 



of the Is — Od shell model using the realistic effective interaction of Wildenthal ||41|| . As 
we have already mentioned, this interaction has been tested over many years. It is known 
to accurately reproduce many nuclear observables for s-d shell nuclei. The Wildenthal 
two-body matrix elements as well as the single particle energies are determined by least 
square fits to experimental data in the region of the periodic table with A = 17 — 39. 
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The eigenstates of the daughter nucleus 27 Na were evaluated in the isospin represen- 
tation. For each spin Jf with T = 5/2 the first states reaching up to E x = 25 MeV in 
excitation energy were calculated. On the other hand, for 27 A1 we evaluated the ground 
state (5/2)+ with T = 1/2. We also constructed all the excited (positive parity) states up 
to 5 MeV in order to check our predictions against experiment. In Fig. 1 we present the 
calculated and measured [42] low-energy spectrum of 27 A1 up to 5 MeV. As can be seen 



from this figure, the spectrum of 27 A1 is well reproduced. In the case of the unstable 27 Na 
isotope the comparison between theory and experiment can not be accomplished due to 
lack of experimental data. 

For the special case of the reaction (|13|) studied in the present work, since M( 27 Al) — 
M( 27 Na) = —10.6 MeV, the momentum transfer at which our matrix elements must be 
computed, is given by 

Pe = {Pe)max ~ E x /c, (p e )max = 94.5 (MeV/c) (47) 



5 Results and discussion 

As we have mentioned, the primary purpose of the present work is the calculation of the 
total fi~ — > e + reaction rate by summing over partial transition strengths. As a first step 
we restricted ourselves in the light neutrino mechanism. As a result we only dealt with 
the Fermi type and the Gamow- Teller type operators discussed in Sect. 2. In other words 



we evaluated the reduced matrix elements [3^j 



M F = (f\\n F \\i{gs)) (48) 

and 

M GT = (f\\Q GT \\i(gs)) (49) 

for the transitions between the initial \i) = (5/2)+ and all the final |/) = |JT) states up 
to 25 MeV. We evaluated the contributions arising from both the real part and imaginary 
part of the radial function F(r) Eq. (|25D which occur in the case of light Majorana 
neutrinos. 

In the calculation of Fermi |Mi?| 2 and Gamow- Teller |Mgt| 2 strengths we found that 
the contribution from the multipolarities L = 2 and L = 4 is, in general, quite small 
compared to that of L = 0. This becomes obvious by glancing at Table 1 where the total 
Fermi and Gamow- Teller strengths with respect to multipolarities L are listed. As can be 
seen from Table 1 the Gamow- Teller strength is almost 9 times greater than the Fermi 
one which means that it dominates the total strength. 

In order to compare the branching ratio originating from the g.s. — > g.s. transition with 
that associated with the transition to all final (5/2) + states, we define, for convenience, 
the ratio 

R gs _ M, 2 



where 



A = '-f- = ^ » ( 5 °) 

S 9 s^f = ( 1 - T^M (51) 

The quantity A gives the portion of the g.s. — > g.s. contribution relative to the total rate 
here computed by the sum over all partial transitions included in our model space. For 
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the g.s. — > g.s. transition p e = 94.5 MeV/c according to Eq. fl4"?|) . Since m e c << p e we 
can consider the approximation p e « -E e /c, which is equivalent to neglecting the electron 
mass (m e ) in the kinematics of the reaction flT3|). 

For the explicit contribution of the excited states to the branching ratio we present in 
Table 2 the sum J2f <3gs-+f f° r eacn set of excited states of given \Jf). The second column 
of Table 2 corresponds to the contribution from the real part of Eq. (pBD while the third 
column represents the total contribution coming from the real and imaginary parts of Eq. 
(p5p. As can be seen from Table 2, the ground state transition exhausts a large portion 
(41%) of the total yT — > e + reaction rate. As can also be seen from Table 2, the main 
contribution to the total rate comes from the (5/2) + multipole states which contribute 
about 98% of the total rate. The rest, 2%, originates mainly from the (3/2) + , (7/2) + 
and (9/2) + states. Another interesting conclusion stemming out of the results of Table 
2 is the fact that, the contribution coming from the imaginary component of the radial 
function F(r) Eq. fl2"5| ) dominates the total branching ratio. In fact the contribution of 
the real part is about 20 times smaller than the corresponding imaginary one. 

In order to have an insight about the origin of this difference we studied the behavior 
as a function of the relative coordinate r of the following three quantities: 

(i) 7 nn /(r) = -R n oRn'o r 2 b Im{F(r)}, 

(ii) S nn ,(r) = RnoRn'o r 2 b Re{F(r)}, 

where n, n indicate the the nodes in the relative coordinate of the two nucleon wave 
function. Im{F(r)} and Re{F(r)} are the imaginary and real parts, respectively, of F(r) 
[see Eq. (p8[) ]. Since there is no interference between the real and imaginary parts we 
find it convenient to define I nn > with the (-) sign. For calculations within the Od — Is shell 
we make the plausible assumption that the main contribution comes from zero angular 
momentum states in the relative motion. 

(hi) D nn /(r) = R n oR n 'o r 2 b, which expresses the nuclear densities and corresponds to 
putting F(r) « 1 (this is the case of 0i//3/3-decay). 

R n i denotes the radial part of the harmonic oscillator wave functions entering the 
matrix elements of Eqs. (^8|) and ( (49]) . Thus, D nn r, I nn > and S nn > are dimensionless 
quantities of which the squares determine the magnitudes of the two-body matrix elements 
involved in our evaluations of the partial transition rates. The results obtained for D nn r(r), 
Inn>{r) and S nn i(r) with n,n' = 0,1, are plotted in Fig. 2 from where we conclude the 
following: 

The variation of I nn >(r) which contains the imaginary part of F(r) is in all cases in 
phase with D nn >(r). The peaks of these quantities appear at about the same value of r for 
each case. On the contrary, the variation of S nn r(r) which contain the real part of F(r), is 
not always in phase with D nn i{r) and I nn i{r). Also the peaks of S nn i(r) do not appear at 
the same positions of r as those of D nn r(r) and I nn >(r). This picture appears in all other 
cases resulting by using the various R n i wave functions entering our calculations, even 
though the absolute strengths of D nn i(r), I nn i{r) and S nn >(r) are much smaller compared 
to those of Fig. 2. 

Obviously, the area bounded by the r-axis and I nn > or S nn i gives a measure of the 
contribution of the imaginary or real part of the /i~ — > e~ operators, respectively, A 
rough estimation of this area gives for the ratio of imaginary to real a value of about 4-5 
which justifies the factor of about 15-25 between imaginary and real part contributions 
to the partial sum evaluation of the total fi~ — > e~ rate. Of course, this difference of the 
imaginary part need not apply when the model is expanded outside our space, i.e. to 
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include lfwo and 2hu excitations. 

At this point we found interesting to study the spreading of the contributions through 
the excitation spectrum of the final nucleus. To this end, in Fig. 3 we plot the distribution 
of Sg S ^f for all (5/2) + states. Similar pictures are obtained for the other multipole states 
(1/2)+, (3/2)+, (7/2)+, (9/2)+, (11/2)+, (13/2)+, even though their contribution is quite 
smaller than that of (5/2) + . The common feature of these distributions is the fact that for 
each multipolarity the main contribution comes from low-lying states and that excitations 
lying above ~10 MeV contribute negligible quantities. In some cases, as e.g. (7/2) + the 
contribution comes from a narrow window of the excitation spectrum of the final nucleus. 

At this stage we should mention that the total rates can also be obtained in the context 
of the closure approximation. We remind that according to the closure approximation the 
contribution of each individual state is effectively taken into account by assuming a mean 
excitation energy E = (Ef) — E gs , and using the completeness relation J2f \ f)(f\ — 1- 
Therefore 

j2\(f\m 2 = m + m 

f 

The matrix element can be written as a sum of two pieces: a two-body term 

and a four-body one, that is 

= (i\(n + n) 2b \t) + (i\(n + n) 4b \i) (52) 

The total rate evaluation thus involves only the g.s. of the initial nucleus. The 4-body 
piece was neglected in the earlier calculations. Thus, following the work of Ref . |19[ , the 



relevant two-body operator for the process studied in the present work takes the form: 
n + n = £M^-(j)T + (z)r + (j)[§ - a(i).a(j)\[% ~ <i).a{j)]A 2 (53) 

i#j J A J A r ij 

The above equation can be written as follows: 

Q+n = Sr.(i)r + (i)r-(j)r + (j)[9 + $ ? - 2(& + l)cr(i).a(j)A 2 (54) 

i#j J A J A r tj 

To proceed further one makes the simplifying assumption that the matrix element is 
dominated by the spin singlet states. Then the isospin operator counts the number of 
proton pairs in the initial nucleus. The radial part is estimated by assuming a uniform 
nuclear two body density, i.e. 

where r c is the hard core radius, assuming for simplicity that the short range correlations 
are described by a simple step function. In any case since r c R , the short range 
correlations can be neglected and to a good approximation the value of the above integral 
is 1/4. We thus find that 



( l \n + n\t) = z(z-i)(£t + 3) 2 j (56) 

J A 4 
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The above matrix element must be multiplied by the kinematical factor [(pe)/{Pe)max\ 2 ~ 
0.8, taking an average excitation energy of about 20 MeV. We thus find: 

« 0.8 x 13 x 12 x (3.7) 2 x 0.25 = 430 (57) 

It is clear, therfore, that the contribution of the excited states found by our state-by-state 
calculation in the present work is much smaller than the predictions found previously by 
employing closure approximation as outlined above. 

The disagreement appeared between closure approximation and the present state-by- 
state calculation can be attributed to the following reasons: 

i) The closure approximation takes into account not only the contribution of Ohoo space 
but also includes excitations E > nhu, n > 1, as well as the continuum spectrum. A 
possible extension of the s — d model space to include two-particle two-hole states with 
the above harmonic oscillator excitations is practically impossible. 

ii) In closure approximation the second term in Eq. (|52|) which includes the four- 
body forces and which is very complicated, was not taking into account in the previous 
calculations. Of course, the obvious question raised is: how important the contribution 
of four-body terms is ? 

iii) From the ordinary /i-capture it is known that the results of the simple closure 
approximation are quite sensitive to the assumed mean excitation energy E. Since the 
spectrum of the final nucleus reached by the operators of — > e + conversion has not yet 
been observed, E may be naively estimated from the spectrum of the final nucleus. 

Anyway putting our nuclear matrix elements and all the other nuclear physics input 
into Eq. flU]) we get 

R pe+ =2Ap\r] eff \ 2 . (58) 
Combining this with the present experimental limit 



R„ e + < 4.4 x lO - " p 
we get 

p\Veff\ 2 < 1-8 x 10~ 12 (59) 

The quantity p|f? e //| 2 depends on the particle model adopted as well as the prevailing 
mechanism for this process and it will not be further discussed in this work. 



6 Summary and Conclusions 

In the present work we have investigated the exotic double-charge exchange neutrinoless 
muon-to-positron conversion in the presence of nuclei, (A,Z){i~ — > e + (A,Z — 2). The 
appropriate operators have been constructed considering a variety of gauge model pre- 
dictions, not only light intermediate Majorana neutrinos. We have chosen to study the 
experimentally interesting nucleus 27 A1, since it is going to be used as a stopping target in 
the Brookhaven experiment. This nucleus is expected to be described rather well within 
the Is — Od shell model, since one can use the well-tested Wildenthal Is — Od interaction. 

As a first step we restricted to the calculation of the rates in the case of light inter- 
mediate neutrinos. 

^From our results on the reaction 27 Al(/j,~, e + ) 27 Na, we can conclude the following: 
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(i) The contribution coming from the Gamow-Teller component of the \i~ —>■ e + 
operator dominates the total rate matrix elements. 

(ii) The contribution of L = multipolarity dominates the total rate. 

(iii) The total strength, resulting by summing over partial transition matrix elements 
included in our model space, is much smaller than that found previously by using closure 
approximation. 

(iv) The contribution arising from the imaginary component of the radial part in the 
pT — ► e + conversion operator is much larger than the one of the real component. 

(v) The portion of the g.s. —>■ g.s. contribution (which is proportional to the matrix 
element |.M gs _> gs | 2 ) into the total rate is 41%. This is good news, since, eventually, the 
experiments will focus on the ground state. 

(vi) By putting our nuclear physics input into Eq. ( ]4"4| ) we obtain the limit p\i] e ff\ 2 — 
1.8 x 10~ 12 . The specific limits on p and rj e ff depend on the particle model assumed and 
the prevailing mechanism for the \i~ — > conversion process. 

Appendix 

According to the gauge models mentioned in Sect. 2 the transition operator Q can be 
given in terms of the following components : 



Vsa = $>_(i)r_(j)e tPer< f(r tJ )W Sa (tj), a = 1,2 
n Aa = Y,T-(*)r-ti)e tPer > f(r t] )W Aa (ij), a = 1,2 

i^3 



where 



Wsx{ij) = 1 
W S 2(ij) = 0i ■ Oj = -\/3 [o-j <g> (Tj 

W S3 (ij) = 3(ai ■ hj)(vj • hj) ~ <*% ■ <*j = VQ [V4KY 2 (hj) <8 

W A1 (ij) = iOi x Oj = (-V2) [a, <g> cxj] 1 
W A2 (ij) = Oi - Oj 

W A3 {ij) = fa - aj^thj x %) = \ l llV^Y 1 ^) ® ^Y 1 {R ij ) 



and VL Ai read 



(Ji (X) a 4 



Oi — Oa 



(60) 
(61) 
(62) 
(63) 

(64) 
(65) 

(66) 

(67) 
(68) 

(69) 



By Applying the usual multipole decomposition procedure the operators f2 Sa , Vt Aa , Vt T 

n Sa = E ^y\Pe) ■ Off )J 5 kL 5 L j (70) 

LA.J 



V Aa = J2 [V^Y A (p e 



o 



{L,S)J 
A a 



LJA 



5 



AL 



(71) 
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Q' T = ]T V^Y A (p e ) ■ 0^ S)J 6 AJ (72) 

LAJ 

n' A2 = £ v^Fy A (p e ) • o^' 5)J 5 AJ (73) 

LAJ 

The operators 0^' 5 ' )J , O^f^ , 0^^ J and ' 5 ' )J are given by the following equations 



IC i<j Z 

[V^Y l (fij) <g> V4^Y c (R ij )] J , ^ = 0, J = L 



(74) 



IC i<j z 



, ,S = 0,J = L 



47rF'(f lj ) <g> v / 47r £ ( J R, i ) <g) (->/3) ® <7j]° 



(75) 



IC i<j Z 



, S=1,J = L, |L±1| (76) 



0£* jJ = E^E^(^)r-(j)/(^)j'K^)jz:(Pe^) 



, S = 1,J = L, \L±1\ 



AttY 1 (fij ) <g> V 4ttF l (Rij ) <g> (<Ti - <Tj) 



(77) 



4 L,5)J = E^E^-(0r-0-)/'(^)iK^)^(Pe%) 



,A = Z, \l±2\,S = 2,A = J (78) 



47rF A (fy) <g> \ / 4:iTY L (R ij ) <g> [a { <g> a 



where 



n '{L,S)J 



E E E ©^L E T^r-umrjM^jcbeRij) 



IC k\ M 



i<j 



AnY^fij) ® V^irY^Rij) 



, 5=1,7 = A 



(79) 



^ICXL - D tc\L A lCA 



(80) 
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u ic\l — 



\l+C+A 



30LA/< 



2 I A 
£ L A 



A 2 L 
2 A 



(20/0|A0) 



(81) 



C\AIC _ a jAIC 



(82) 



^ = 3v^(-l) i+A+1 ^V / S^A(10/0|A0)(£010| /U 0)( \ \ \ 





/ A 1 




£ 




u 


A; A J 


f; 


L 





A 


1 


fc 


1 


L 


1 


A 


1 


L 


1 


A 






(83) 



with 



l£ 



-^(/0£0|A0)(1 + (-iy)t 



and 



B, 



ICA 



l£J 



\ 3A 



,J+i 



(Z0£0|A0)(1- (-iy)i 



l\J+C 



(84) 



(85) 



(a = 2a + 1). 

We are now going to discuss the operators entering the leptonic R-L interference and 
in some SUSY mechanisms. Now the relevant operators may have a time component, 
which is small and, in any case, except for their radial part is the same with the F and 
GT discussed above. They also have a space component, which is proportional to f\j. 
They are vectors, which yield a scalar, when combined with the leptonic current. They 
are of the form 



' r ij, 



{<7i ■ r i3 )oj - {pi ■ <Jj)rij and %{ai - Oj) x 



The above operators can also be written as 

°A a 3 ■ f'ij) + (°i ' - ( a i ■ VjYij = ^'(0) + c/(2) 



and 
where 



i{<Ji - Gj) X fij = J (I) 

J(k) = a{k) [VA^Y 1 ^) <g> T k (spin)}\ k = 0, 1, 2 
1 



a(0) 



a(2) 



3a/3 
2^ 



, T°(spin) = Oi ■ Oj 



, T 2 (spin) = [(Tj (g) <jj 



The above operators are accompanied by the lepton outgoing waves 

Oij = exp(ip e ■ Rij)[exp(ip e ■ r^)/2 + (-l) fc+1 exp(-ip e ■ r ij )/2] 



(86) 

(87) 

(88) 
(89) 

(90) 

(91) 

(92) 
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The phase of the second term guarantees that the combined operator is overall symmetric 
under the exchange of the particles % and j. The last operator can be brought into the 
form 



On = ]T [(3{lCkA)j l {p e r ij /2)j c {p e R ij )[[^Y\^ (93) 

ICXA 

where 

(5{lCkA) = [1 + (-l)' +fe+1 ] t l+c (-l) l+c [(2l + 1)(2£ + 1]^ 2 (IOJCO\\0) (94) 
Combining the above factors we obtain 

Q'(k) = J2 [V^Y A (p e ) <g> (L ' 5)J (A;)] 1 (95) 

LJA 

where 

0^ J (k) = a(k) £ r_(i)r_(j)/(r«) £ W^Hh \ C, L, k, J, A) j,(p e r y /2) 

j £ (p e %) [[v^Fy A (r^-) ® v / 4^F £ (%)] L ® T fc (spin)] J (96) 

with 

-y(!,\,£,L,h,J,A) = [(3(2L + 1)(2/ + 1)(2J+1)] 1 / 2 (97) 



(lOZOjAO) 



1 k 1 1 J 1 I X 
J A L M C L A 
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Multipole Fermi Contr. Gamow-Teller Contr. 



L 


= 


8.291 


77.272 


L 


= 2 


0.076 


0.426 


L 


= 4 


1.553 xlO" 4 


5.515 xl0~ 4 



Table 1: Fermi and Gamow-Teller transition strengths for various multipolarities con- 
tributing to the partial rate of all the (5/2) + states. Both the real and imaginary parts 
of Eq. (|25|) are included. 



1/) 


Real 


Total 


A (%) 


gs -> gs 


1.38 


31.59 


41.44 


1/2 


3.63xl0- 3 


0.044 


0.05 


3/2 


6.09xl0- 2 


0.55 


0.72 


5/2 


2.24 


43.52 


57.08 


7/2 


3.39xl0~ 2 


0.29 


0.37 


9/2 


3.37xl0- 2 


0.24 


0.31 


11/2 








13/2 








Total 


3.75 


76.23 


100 


Closure Approx. 




430 


100 



Table 2: Individual contributions of J2f Sgs^f arising from the real and imaginary part of 
the integral F(r) of Eq. (j25]) . As we can see, the process is dominated by the contribution 
of the imaginary part. 
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FIGURE CAPTIONS 



FIGURE 1 The calculated (right) and measured (left) |42| energy spectrum for the 
lowest positive parity states of 27 Al. 

FIGURE 2. Variation of the quantities D nn r, I nn i and S nn i as functions of the relative 
coordinate r (see Sect. 5) assuming that the relative zero angular momentum states 
dominate. Since the relative motion is indicated in the plots one sees some contribution 
beyond the nuclear radius. In this figure only the most prominent cases for n, n' = 0, 1 
are shown. The behavior of D nn >(r), I nn '(r) and S nn i{r) in the other cases needed for our 
calculations is similar. 

FIGURE 3. Distribution of the transition strength S gs ^f for all J n = (5/2) + states. 
Both the real and imaginary parts of Eq. (|25|) are included. 
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